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We review the current status of theoretical study of non-leptonic two body B 
decays. There are two independent directions for this purpose. One is the so 
called QCD factorization approach (or BBNS approach), which is based on naive 
factorization approach. The other one is named perturbativc QCD approach. We 
list the different ideas and applications of the two method, and make a comparison 
between the two. 

1 Introduction 

The current running B factories at KEK and SLAC arouse many interests of 
theoretical studies of B decays. The mechanism of CP violation and signal of 
new physics are among the most important issues in B physics. Most of B topic 
involve the study of non-leptonic B decays, such as CKM angle measurements 
and rare decays etc. Quantum Chromo dynamic is one of the successful theory 
of particle physics, however, the non-pcrturbative behavior of QCD is still 
the unsolved problem. Unlike perturbativc QCD, the hadronization of non- 
lcptonic B decays can not be done from the first principal. Thus it is model 
dependent. The most successful approach recently is the factorization approach 
), which can explain many of the decay channels by very few parameters 



Although the predictions of branching ratios agree well with experiments 
in most cases, there are still some theoretical points unclear in FA. First, it 
relies strongly on the form factors, which cannot be calculated by FA itself. 
Secondly, the generalized FA shows_that the non-factorizable contributions are 
important in a group of channels aa. The reason of this large non-factorizable 
contribution needs more theoretical studies. Thirdly, the strong phase, which 
is important for thx; CP violation prediction, is quite sensitive to the internal 
gluon momentum □. This gluon momentum is the sum of momenta of two 
quarks, which go into two different mesons. It is difficult to define exactly in 
the FA approach. To improve the theoretical predictions of the nondeptonic 
B decays, we try to improve the factorization approach, and explain the size 
of the non-factorizable contributions in a new way. 

There are mainly two directions toward this improvement: One is the so 
called QCD factorization or BBNS til, which is based on the naive factorization 
approach. In this approach, the non-factorizable contribution can be directly 
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calculated if it is not dominant over the factorizable one. The other is the 
pertuxbative QCD (PQCD) approach which is based on Brodsky and Lepage's 
idea U. In this direction, people try to calculate the form factors and also 
the non-factorizable and annihilation type contributions in a systematic way. 
In the next section, we will first explain the idea of naive and generalized 
factorization approach, which has been studied for years. In section 3, we will 
introduce the idea of QCD factorization. And in section 4, PQCD approach is 
applied. Finally, in section 5, we will give a brief comparison between the two 
approaches and summarize at the end. 



2 Generalized Factorization Approach 

The calculation of non-leptonic decays involves the short-distance Lagrangian 
and the calculation of hadronic matrix elements which are model dependent. 
The short-distance QCD corrected Lagrangian is calculated to next-to-leading 
order. The popular method to calculate the hadronic matrix elements is using 
the factorization method where the matrix element is expressed as a product 
of two factors (hih2\H e f f\B) = (hi\ Ji\B) (Ji2\J2\0) ■ The first factor is propor- 
tional to the B — ► hi form factor, while the second one is proportional to the 
decay constant of /12 meson. 

First, let us discuss the short distance part. The non-leptonic decays of B 
mesons are induced by the weak interaction in the quark level, which gives the 
effective four-quark operator. Together with the QCD corrections, the effective 
Hamiltonian for the charmless non-leptonic B decays is 
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where q = d,s and V q > q denotes the CKM factors. The operators 0\,02 
are tree level current operators. The operators O3,... ,Oq are QCD penguin 
operators. O7, ...,Oiq arise from electroweak penguin diagrams, which are 
suppressed by a/a s . Only O9 has a sizable value whose major contribution 
arises from the Z penguin. 

For example, let us consider B meson decays to two pseudoscalar mesons. 
In the factorization approach, using the effective Hamiltonian, we write the 
required matrix element in its factorized form 

(P^lHefflB) = l ^£v qb V; q ,a l f P2 (m 2 B - m\)F^{ml). (2) 

The dynamical details are coded in the quantities Oj, which we define as 
Oi = Cf f + Cf f /N c , where {i,j} is any of the pairs {1,2}, {3,4}, {5,6}, 
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{7,8} or {9, 10}. In practice, N c is treated as a phenomenological parameter 
to include the non-factorized color-octet contributions. 

If the amplitude is dominated by the tree amplitude, the BSW-classification 
can be appliectl. Class-I decays are color favored, whose matrix elements are 
proportional to a\ = C\/N c + C2. Class-II decays are color suppressed, whose 
matrix elements are proportional to a-2 — C\ + C2/N c . If N c = 3 as in QCD, 
then ci2 becomes a very small number, making the branching ratios of this 
class of decays also small. However from experiments we know that in this 
category of decays, the branching ratios are not so small. Therefore the non- 
factorizable contributions are found to be very important in order to explain 
the large branching ratios. Class-Ill decays are proportional to ai + rci2- This 
class of decays will determine the relative sign between a% and ci2- The QCD 
and electroweak penguins are also present in the charmless decays oiB mesons. 
For the penguin-dominant decays, we introduce two more classes B: Class-IV 
decays involve one or more of the dominant penguin coefficients 04, ciq and ag. 
Class- V decays are decays with strong 7V c -dependent coefficients 03, as, 0,7 and 
aw 

Class-I and Class-IV decays have relatively large branching ratios of the 
order of 10~ 5 and stable against variation of N c . Most of the measured decay 
channels by experiments are belong to these classes. There is a good agreement 
between experiments and theory, which indicate the success of FA. Class-Ill 
decays are mostly stable, except for some B — > PV decays. Class-II and Class- 
V decays are rather unstable against variation of N c . Some measured decay 
channels of these classes indicate that we need an effective N c ~ 2 rather 
than N c — 3 to explain the experimental data. This implies that large non- 
factorizable contribution exist in these decays. Many of them may receive 
significant contribution from the annihilation diagrams and/or soft final state 
interactions. 

In FA, the strong phase comes only from the so called BSS mechanismi. 
The inner charm quark loop produce a strong phase, when the charm quarJc 
on mass shell. This strong phase depends on the inner gluon line strongly C. 
Therefore, in FA, the predicted CP violation of B decays are in a quite large 
range. 

In the study of factorization approach, one finds that the non-factorizablc 
contributions are important in some of the non-leptonic B decays, which can 
not be explained in the factorization approach itself. Therefore theoretical 
study of this effect is required. 
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3 QCD Factorization Approach 

Recently, Beneke, Buchalla, Neubert, and Sachrajda (BBNS) proposed a for- 
malism for two-body charmless B meson decays In this approach, they 
expand the hadronic matrix element by the heavy b quark mass 

(7T7r\Q\B) = (7r|ii|fl)<7r|j 2 |0) [l + ^r n < + G(A QC D/m b )] , (3) 

where Q is a local operator in the effective Hamiltonian and ji^ are bilinear 
quark currents. By neglecting the power corrections in Aqcd, one need only 
calculate the order a s corrections including the vertex corrections for the four 
quark operators and the non-factorizable diagrams. These diagrams are shown 
in Fig.l. The first 6 diagrams have already been included in the generalized 
FA approach as next-to-leading order QCD corrections to local four quark 
operators. What new are the last two non-factorizable diagrams, which has a 
hard gluon line connecting the four quark operator and the spectator quark. 

They claimed that factorizable contributions, for example, the form factor 
F B * in the B — > 7T7T decays, are not calculable in PQCD, but nonfactorizable 
contributions are in the heavy quark limit. Hence, the former are treated in 
the same way as FA, and expressed as products of Wilson coefficients and 
F Bn . The latter, calculated as in the PQCD approach, are written as the 
convolutions of hard amplitudes with three (B,tt,tt) meson wave functions. 
Annihilation diagrams are neglected as in FA. Hence, this formalism can be 
regarded as a mixture of the FA and PQCD approaches. Values of form factors 
at maximal recoil q 2 = and nonperturbative meson wave functions are all 
treated as input parameters. It is easy to see from eq.(||), that this equation is 
only applicable for those color enhanced decay modes, where the factorizable 
contribution dominates the final results. While for the color suppressed modes, 
where the non-factorizable contributions are not small, the expansion of eq.([|) 
is not right, since the large non-factorizable contribution is grouped into the 
second term of eq. (||) . 

When extending the BBNS formalism to the B — > D^ir decays, difficulty 
occurs in the calculation of nonfactorizable amplitudes. The decay channel 
of B° — > ir + D~ is calculable, because the end-point singularities in the non- 
factorizable diagrams (last diagrams in Fig.l) cancel each other. However, 
such a soft cancellation does not occur in the calculation of B + — > ir + D° and 
B° — > ir°D° decays. The D meson is heavy on the upper side of the diagram, in 
which the light quark and the heavy c quark do not move collinearly. Therefore, 
soft gluons can resolve their color structure, and interact with them. That is, 
not all the nonfactorizable amplitudes in the BBNS formalism are calculable. 

There are many calculations for various decay channels in this approach. 
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Figure 1: Figures calculated in the QCD factorization approach. 



Slroup of people calculate B meson decays to two light pseudoscalar mesons 
, B meson decays to one pseudoscalar and one vector meson □ and B meson 
decays to final states with r\ or rf E3 in the QCD factorization approach. The 
numerical results show that the theory and experiments agree well for those 
class I and IV decays, which are color enhanced and dominated by the factor- 
izable contribution. This also agrees with the FA result, since the dominant 
part in eq.(||) is the same as the FA. The success of QCD factorization is that 
one can calculate the sub-leading 0{a s ) non-factorizable contribution (second 
term in eq.(Q)) using perturbative QCD. While in the FA, one has to input a 
free parameter N^* to accommodate the non-factorizable contribution. 

Cheng and Yang did the calculation of B meson decay to two vector mesons 
til, including B — > <pK* . In the calculation of B — > <f>K decay, they found that 
the annihilation contributions are important, but the relative strong phase,can 
not be predicted exactly, due to the sensitivity of the cut-off introduced E3. In 
fact, they found that a real annihilation contribution is required to enhance 
the B — ► <f)K decay branching ratio. In the calculation of B — > J/tpK decay, 
they found that the leading twist contribution is too small to accommodate 
the experimental dataL3. Problem remains for color suppressed decay modes 
in QCD factorization approach. 

In the calculations above, those people found that there exist endpoint 
divergence in the annihilation diagram calculations of QCD factorization 
Logarithm divergence occurred at twist 2 contribution, and linear divergence 
exists in twist 3 contribution. If not symmetric wave function, like meson, 
there is also soft divergence in the non-factorizable diagrams. It is very difficult 
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to treat these singularity in the BBNS approach. A cut-off is introduced to 
regulate the divergence, thus makes the QCD factorization approach prediction 
parameter dependent, especially for the strong phase. Recently an effort is 
made to introduce the fey dependence of the wave functions, and Siidakov 
form factors in the BBNS approach in order to remove the singularitiesE£l. This 
makes the BBNS approach to go toward the direction of PQCD approach. 

As for the strong phase in BBNS, like in FA, it may come from the BSS 
mechanism. Here the momentum of the inner gluon is well defined. However, 
it predicts too small strong phase, because of the small gluon momentum. 
There is also another source of strong phase from the annihilation diagrams, 
but strongly depends on the cut-off parameter. The strong phase in QCD 
factorization can be almost arbitrary large. 

Finally, the QCD factorization approach is at least one step forward from 
Naive Factorization approach. It gives systematic prediction of sub-leading 
non-factorizable contribution for the class I and class IV decays, which are 
dominated by the factorizable contribution. Big problem is the endpoint sin- 
gularity, but may be solved with Sudakov from factors like PQCD approach. 
The input parameters in QCD factorization approach are form factors, wave 
functions etc. 

4 Perturbative QCD Approach 

In this section, we will introduce the idea of PQCD approachQ. The three 
scale PQCD factorization theorem has been developed for non-lepie-nic heavy 
meson decay&HI, based on the formalism by Brodsky and Lepage B, and Botts 
and StermantS. In the non-leptonic two body B decays, the B meson is heavy, 
sitting at rest. It decays into two light mesons with large momenta. Therefore 
the light mesons are moving very fast in the rest frame of B meson. In this case, 
the short distance hard process dominates the decay amplitude. The reasons 
can be ordered as: first, because there are not many resonance near the energy 
region of B mass, so it is reasonable to assume that final state interaction is 
not important in two-body B decays. Second, With the final light mesons 
moving very fast, there must be a hard gluon to kick the light spectator quark 
(almost at rest) in the B meson to form a fast moving light meson. So the 
dominant diagram in this theoretical picture is that one hard gluon from the 
spectator quark connecting with the other quarks in the four quark operator 
of the weak interaction. Unlike the usual factorization approach, the hard 
part of the PQCD approach consists of six quarks rather than four. We thus 
call it six-quark operators or six-quark effective theory. There are also soft 
(soft and collinear) gluon exchanges between quarks. Summing over those 
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leading soft contributions gives a Sudakov form factor, which suppresses the 
soft contribution to be dominant. Therefore, it makes the PQCD reliable in 
calculating the non-leptonic decays. With the Sudakov resummation, we can 
include the leading double logarithms for all loop diagrams, in association with 
the soft contribution. 

There are three different scales in the B meson non-leptonic decay. The 
QCD corrections to the form-quark operators are usually summed by the renor- 
malization group equation cB. This has already been done to the leading log- 
arithm and next-to-leading order for years. Since the 6 quark decay scale rrif, 
is much smaller than the electroweak scale mw, the QCD corrections are non- 
negligible. The third scale 1/6 involved in the B meson exclusive decays is 
usually called the factorization scale, with 6 the conjugate variable of parton 
transverse momenta. The dynamics below 1/6 scale is regarded as being com- 
pletely non-perturbative, and can be parameterized into meson wave functions. 
The meson wave functions are not calculable in PQCD. But they are univer- 
sal, channel independent. We can determine them from experiments, and it 
is constrained by QCD sum rules and Lattice QCD calculations. Above the 
scale 1/6, the physics is channel dependent. We can use perturbation theory 
to calculate channel by channel. 

Besides the hard gluon exchange with the spectator quark, the soft gluon 
exchanges between quark lines give out the double logarithms ln 2 (P6) from 
the overlap of collinear and soft divergence, P being the dominant light-cone 
component of a meson momentum. The resummation of these double loga- 
rithms leads to a Sudakov form factor exp[— s(P, 6)], which suppresses the long 
distance contributions in the large 6 region, and vanishes as 6 > 1 /Aqcd ■ This 
form factor is given to sum the leading order soft gluon exchanges between the 
hard part and the wave functions of mesons. So this term includes the dou- 
ble infrared logarithms. It is shown in refE3 that e~ s falls off quickly in the 
large 6, or long-distance, region, giving so-called Sudakov suppression. This 
Sudakov factor practically makes PQCD approach-applicable. For the detailed 
derivation of the Sudakov form factors, see refEjO. 

With all the large logarithms resummed, the remaining finite contributions 
are absorbed into a perturbative b quark decay subamplitudc H{t). Therefore 
the three scale factorization formula is given by the typical expression, 



C(t) x H(t) x $(x) x exp 



-s(P,b)-2 f ^ 9 (aM) 
Ji/b M 



(4) 



where C(t) are the corresponding Wilson coefficients, Q(x) are the meson wave 
functions and the variable t denotes the largest mass scale of hard process 
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Figure 2: Diagrams for _B — + P\P^ decay in perturbative QCD approach. The factorizable 
diagrams (a),(b), non-factorizable (c), (d), factorizable annihilation diagrams (e),(f) and 
non- factorizable annihilation diagrams (g),(h). 



H, that is, six-quark effective theory. The quark anomalous dimension 7 g = 
—a s /n describes the evolution from scale t to l/b. Since logarithm corrections 
have been summed by renormalization group equations, the above factorization 
formula does not depend on the renormalization scale fj, explicitly. 

The 7r meson is treated as a light-light system. At the B meson rest frame, 
pion is moving very fast. We define the momentum of the pion which contain 
the spectator light quark as P 2 = (tub / V%) (1) 0, Ot) ■ The light spectator quark 
moving with the pion (with momentum P 2 ), has a momentum (fc^,0, tczr)- If 
we define the momentum fraction as x-i = fc^/P^, then the wave function of 
pion can be written as 



1 



= 75 [^7r07r(x2,k 2 T) + to o 0£ (x 2 , k 2T ) + TO O 0£ (x2 , kar)] , (5) 



where 4>n{x2i k 2 r) is twist- 2 wave function and (x 2 , k 2 ^) and are twist-3 
wave functions. The toq m ec l-(|l) is given as 



mo 



(G) 



It is not the pion mass. Since this too is a scale characterizing the Chiral 
symmetry breaking, it is estimated around 1 ~ 2 GeV using the quark masses 
predicted from lattice simulations, one may guess contributions of Too term 
cannot be neglected because of too tub- In the previous calculation of 
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B — > imtA and _B — > 7r_K~ decays E3 we do not include the last term in eq.(^|) 
for the pion wave function. However, by using a phenomenology twist 3 wave 
function for we get the right result for those branching ratios. The reason 
is that this choice of twist 3 wave function accommodate the full twist 
3 contribution effectively. This phenomenological model also accommodate 
effectively the threshold resummation effect discussed below. 

In the PQCD approach, we can calculate not only the factorizable dia- 
grams (Fig. 2(a), (b)) and non-factorizable diagrams (Fig. 2(c), (d)) contribution 
but also the annihilation type diagrams (Fig.2(e,f,g,h,)). Unlike the QCD 
factorization approach, there is no logarithm or linear divergence in these cal- 
culations due to the reason of Sudakov suppression with kx resummation and 
the threshold resummation discussed above. 

As shown above, in the PQCD approach, we keep the kx dependence of the 
wave function. In fact, the approximation of neglecting the transverse momen- 
tum can only be done at the non-endpoint region, since kr <C k + is qualified 
at that region. At the endpoint, k + — > 0, kx is not small any longer, neglecting 
kx is a very bad approximation. By, keeping the kx dependence, there is no 
endpoint divergence as occurred in the QCD factorization approach, while the 
numerical result does not change at other region. Furthermore, the Sudakov 
form factors suppress the endpoint region of the wave functions. Recently an- 
other type of resummation has been observed. The loop, correction to the weak 
decay vertex produces the double logarithms a s In 2 xj£^- Using the wave func- 
tions from light-cone sum rules, at the endpoint region, these large logarithms 
are important, they must be resummed. The threshold resummation for the 
jet function results in Sudakov suppression, which decreases the contribution 
of endpoint region of wave functions. Therefore, the main contributions to the 
decay amplitude in PQCD approach comes not from the endpoint region. The 
perturbative QCD is applied safely. 

By including kx to regulate the divergence, large logarithmic corrections 
a s In kx appear, and Sudakov resummation is demanded. With the resultant 
Sudakov suppression, we have explicitly shown that almost 100% of the full 
contribution to the B — ► n transition form factor arises from the region with 
the coupling constant a s /ir < 0.3 ca. It indicates that dynamics Jjrom hard 
gluon exchanges indeed dominate in the PQCD calculation. In til Sudakov 
resummation is irrelevant, since all QCD dynamics has been parameterized 
into models of form factors. 

We emphasize that nonfactorizable and annihilation diagrams are indeed 
subleading in the PQCD formalism as Mb "fh? 6 ^ This can be easily observed 
from the hard functions in appendices of refEa-El When Mb increases, the B 
meson wave function enhances contributions to factorizable diagrams. How- 
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ever, annihilation amplitudes, being independent of B meson wave function, 
are relatively insensitive to the variation of Mb • Hence, factorizable contribu- 
tions become dominant and annihilation contributions are subleading in the 
Mb — > oo limit El Although the non-factorizable and annihilation diagrams 
are subleading for the branching ratio in color enhanced decays, they provide 
the main source of strong phase, by inner quark or gluon on mass shell. The 
BSS mechanism strong phase is negligible in the PQCD approach. In fact, the 
factorizable annihilation diagrams are Chirally enhanced. They are not negli- 
gible in PQCD approach EjO. In the decays of B meson to two light mesons, 
we collect terms up to chiralltt enhanced terms (O(mo/ms)), but still drop the 
terms suppressed by A/msEa. 

The main input parameters in PQCD are the meson wave functions. It is 
not a surprise that the final results are sensitive to the meson wave functions. 
Fortunately, there are many channels involve the same meson, and the meson 
wave functions should be process independent. In all the calculations of PQCD 
approach, we follow the rule, and we find that they can explain most of the 
measured branching ratios of B decays. FOjL.example: B — > 7i"ZL_decays E3, 
B — > K% decaysc3, B — > up, B — > 7TLa-decayscJ, B — > ATlfj-cLecays ell, the form 
tor calculations of B — » n, B — > pEa, B — > Krfl) decayscj, B — » K(f> decays 



5 Summary 

For a comparison, in the QCD factorization: Form factors are input parame- 
ters, which are claimed to be dominated by soft contribution and not calcu- 
lable. The endpoint singularity is a crucial point in this approach, a cut off 
is needed to regulate the divergence. The QCD factorization follows FA, in 
which it has been assumed that factorizable contributions, being the dominant 
contribution. And all other contributions such as non-factorizable and annihi- 
lation diagrams, being a s corrections, are sub-leading. In the limit of a s — ► 0, 
it goes back to FA. Therefore the current BBNS approach can not be applied 
to the non-factorizable dominant process such as B° — ► D°tt°, and also those 
annihilation diagram dominant processes. 

In the PQCD approach, the form factors are calculable, which are dom- 
inant by short distance contribution. By including the hp dependence and 
Sudakov suppression, there is no endpoint divergence. In the PQCD formal- 
ism non-factorizable amplitudes are of the same order as factorizable ones in 
powers of 1/Mb, which are both 0(1/(MbAqcd))- The smaller magnitude 
of nonfactorizable amplitudes in color enhanced decays are due to the can- 
cellation of the two non-factorizable diagrams. From the viewpoint of power 
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counting, they are of the same order. In case of B — ► Dir decays, the cancella- 
tion is absent. The power counting changes, so that we can also calculate the 
non-factorizable contribution dominant processes. 

For a conclusion, non-Leptonic B decays are important in B physics. The 
theory of non-leptonic B decays is a challenging work. There are still problems 
to be solved. The QCD factorization and PQCD approach are going toward 
the same object, although there are still some differences. They will be tested 
by experiments soon. 
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